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We examine the possibility that the superconductivity in the newly discovered FeAs materials 
may be caused by the Coulomb interaction between d-electrons of the iron atoms. We find that 
when the Hund's rule ferromagnetic interaction is strong enough, the leading pairing instability is 
in spin-triplet p-wave channel in the weak coupling limit. The resulting superconducting gap has 
nodal points on the 2D Fermi surfaces. The k dependent hybridization of several orbitals around a 
Fermi pocket is the key for the appearance of the spin-triplet p-wave pairing. 



I. INTRODUCTION 

Recently, a new class of superconductors - iron based 
superconductors - was discovered. The superconduct- 
ing transition temperature can be as high as 52K.® The 
undoped samples (for example LaOFeAs) appear to have 
a spin ordered phase below 150K.^'^'^*' The electron 
doped^ LaO 1-2; Fa; FeAs and hole doped^ LaOi_a;Sra;FeAs 
samples are superconducting with Tf. ^ 25K. The mag- 
netic field dependence of the specific heat in the electron 
doped material suggests the presence of gapless nodal 
lines on the Fermi surface.'^ 

It appears that the electron-phonon interaction is 
not strong enough to give rise to such high transition 
temperature.^^ In this paper, we will examine the possi- 
bility that Coulomb interaction between c?-electrons on Fe 
drives the superconductivity in the electron/hole doped 
samples and the spin-order in the undoped samples. In 
this case, we find that a p- wave spin-triplet superconduct- 
ing order with gapless nodal lines is the most likely su- 
perconducting order in the weak coupling limit. Naively 
short range repulsion does not have the requisite k de- 
pendence to drive p-wave pairing. It turn out that the 
fc dependent hybridization of several orbitals around a 
Fermi pocket makes the p-wave spin-triplet pairing pos- 
sible. The same model is also shown to have spin-density 
wave (SDW) order for undoped samples. 



II. A THREE-ORBITAL TIGHT-BINDING 
MODEL 

First, let us examine the Fermi surfaces of the iron- 
based superconductor. For concreteness, we will consider 
the LaOFeAs sample. We will assume that the properties 
of the sample are mainly determined by the Fe-As planes. 
The Fe atoms in a Fe-As plane form a 2D square lattice 
(see Fig. 1). Due to the buckling of the As atoms, the 
real unit cell contains two Fe atoms. The real unit cell is 
also a square (see Fig. 1). According to band structure 
calculattions, the Fermi surfaces of the Fe rf-bands are 
formed by two hole pockets at the F point and two elec- 
tron pocket at the M point (see Fig. 2a). ^^"^^ All those 
pockets, mainly formed by the dxz, dyz and dxy orbitals 
of Fe,^^ have similar size, shape, and Fermi velocity. 




FIG. 1: (Color online) (a) The Fe-As plane and the d-orbitals 
of Fe. The filled dots are Fe atoms and the empty dots are As 
atoms. The plus and minus signs in the empty dots indicate if 
the As atom is above or below the Fe plane. The large square 
is the 2D unit cell. The dashed square is the reduced unit cell 
which contains only one Fe. The d^z, dyz and d^cy orbitals of 
the Fe atoms are described by the red, blue and green curves 
respectively, (b) The y hopping between the d^z or dyz with 
d:cy orbitals. 

To understand the mixing of the orbitals near those 
Fermi pockets, we follow Ref. 13,14 to unfold the band 
structure to the extended Brillouin zone (BZ). Ordinar- 
ily, such an unfolding of the band structure extends the 
BZ superficially and there are certain ambiguities in as- 
signing the location of each band. We emphasize that 
this is not the case here. The band structure of the 
extended BZ is uniquely defined due to an additional 
symmetry, i.e. the Fe-As plane is invariant under PzT^ 
and PzTy, where (Ty) is the translation in the x (y) 
direction by the Fe-Fe distance and Pz is the reflection 
z —z (see Fig. 1). Thus, if we combine the translation 
and the reflection P^, then the electron hopping Hamil- 
tonian H has a symmetry described by a reduced unit 
cell with only one Fe per unit cell (see Fig. 1). Since 
[PzTx,H] = [PzTy,H] = [PzTx,P,Ty] = 0, we can use 
the eigenvalues of PzT^ and PzTy to label the single-body 
energy eigenstatcs: 

P,r,|fc) = e'^==|fc), PzTy\k) ^ e''^y\k), (1) 

where fc — {k^, ky) plays a role of crystal momentum. We 
will use such a pseudo crystal momentum to label states 
in an energy band. It is the pseudo crystal momenta k 
that form the extended Brillouin zone (see Fig. 2b) which 
corresponds to the reduced unit cell with only one Fe. 
Let us work out an explicit example by considering a 
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FIG. 2: (Color online) (a) The Fermi surfaces of the Fe ri- 
bands. The plus sign marks the hole pockets and the minus 
sign marks the electron pockets. The square is the Brillouin 
zone, (b) The extended Brillouin zone (dashed square) for the 
reduce unit cell and the positions of the Fermi pockets in the 
extended Brillouin zone. The F point has {kx,ky) = (0,0). 
The folding of the extended Brillouin zone in (b) produces 
(a). After folding, the (f , M) and {X,Y) in the extended 
Brillouin zone map into the F and M in the original Brillouin 
zone, respectively. The yellow shading marks the region where 
the p-wave pairing order parameter may have the same sign. 



tight-binding model involving the d^z-, dy^ and dxy or- 
bitals. First, consider the hopping terms that do not mix 
the orbitals: 
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FIG. 3: (Color online) (a) The zero energy contour of e^zik) 
(red) and eyzik) (dashed blue). The ± are signs of e^zik) 
and €yz{k) in the region. The F point has fc = (b) The 
hybridization of the d^z and dyz bands. The curves are the 
zero energy contours of the hybridized bands, (c) The zero 
energy contour of e^yik). (d) The solid curves are the Fermi 
surfaces of the three-band tight-binding model as a result of 
the hybridization of (b) and (c). 
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where i, j label the positions of the Fe. Since the dxz, 
dyz and d^y orbitals are eigenstates of Pz reflection, the 
PzTx and PzTy symmetries require that t^J, t^j , and t^j 
only depend on i ~ j . The mixing term between the dxz 
and dyz orbitals is given by 
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be the operators with pseudo crystal momentum fe, where 
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where t^^'^^ also depends only on i — j as required by 
the PzTx and PzTy symmetries. On the other hand, the 
mixing term between the d^z and d^y orbitals as well as 
that between the dyz and d^y orbitals are given by 
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where t^^'^^ , and only depend on i — j. We note 

that the dxz and dxy orbitals have opposite eigenvalues 
±1 under Pz- The PzTx and PzTy symmetries require the 
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presence of the factor (—)*==+*!/. Thus, c^^ and c"^ with 
conventional crystal momentum fe - 
with crystal momentum k where Q = (tt, tt) 

However, in the pseudo crystal momentum k space, 
only operators with the same pseudo crystal momentum 



We see that the pseudo crystal momentum k and the 
conventional crystal momentum k are related by fe = fe 



for the dxy orbital and fe 
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k for the dxz and dyz 



orbitals. The total hopping Hamiltonian H ~ Hi 
H3 can be written as H = ^\^Ik^k '^ith 
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It is worth noting that the above Hamiltonian and the 
resulting energy bands are defined on the extended Bril- 
louin zone (labeled by k) of the reduced unit cell (see 
Fig. 2b). 

The nearest neighbor admixture between dxz and dyz 
vanishes by symmetry. Keeping the next nearest neigh- 
bor term only, we have exz,yz = ~^'t'xz.yz['^'^^{kx + ky) — 
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FIG. 4; (a) The energy bands near Y and F. The bands near 
Y come from the d:cz and d^j/ orbitals. The thickness of the 
curve indicates the weight in the d^z orbital. The bands near 
F came from the d^z and dyz orbitals. (b) The contour plot 
of iukVu near Y . The dashed loop is the Fermi surface. 



cos(kx — fcj/)]- Note that this vanishes at F, M, X and Y 
points (see Fig. 3), but is maximal midway between \ 
and Y . Next, consider the nearest neighbor admixture 
between d^z or dy^ with dj.y along the y direction (see Fig. 
lb) . These orbitals can admix only because of the asym- 
metry introduced by the As ions, which has xz symme- 
try. This implies that the dxz-dxy overlap integral is odd 
under x — )■ —x and vanishes. Only the dyz-dxy matrix el- 
ement survives. Also note that under the 180° rotation in 
the x-y plane about the i site, {cl^yYc^^ ~{(^i-yV^i^ ■ 
Thus only the combination {cl^y^c^^ — (cf^y)^c^^ that 
preserves such a symmetry can appear in the hopping 
Hamiltonian which has a form 

i 



k = X and k = Y, and disperses rapid upwards toward F 
but rather flat towards M. The Fermi surface is sketched 
in Fig. 3c. Now we turn on exz,xy and Cyz^xy This gives 
the electron pockets centered at X and Y which become 
the two electron pockets at M after the folding. Note 
that because exz.yz = at X and Y , the band at Y (or 
X) is purely dxz + dxy (or dyz -f dxy). Furthermore, near 
Y, exz,xy ~ sinkx- Thus along Y-T, the hybridization is 
zero and the dxz and dxy bands cross as shown in Fig. 4a. 
Along Y-M, the upper band is mostly dxz while the dxz 
amplitude in the lower band increase linearly with the 
distance from Y . These features are in agreement with 
band calculations. ^^'^^ and are strong evidences that our 
assignment of the dxz band is correct. 

The final result is an ellipse shaped pocket, with the 
long axis of the ellipse pointing in the y direction at Y 
(towards the hole pockets at F). The Fermi surface cross- 
ing along the long direction is purely dxy and the short 
direction is mainly dxz (see Fig. 3d). 

The pocket formed between the dashed lines in Fig. 3d 
are eliminated for sufficiently strong hybridization. How- 
ever, within the 3-orbital model, it is impossible remove 
the Fermi surface surrounding M (the gold solid loop in 
Fig. 3d), because the hybridization elements are all zero 
at M. We need a 4th band which crosses and hybridize 
with the dxz + dyz band to eliminate the unwanted Fermi 
surface. Apart from this, the 4th band plays no role as 
far as the remaining two electron and two hole pockets 
are concerned. It is in this sense that we maintain that 
the 3-orbital model gives an adequate description of the 
low energy Hamiltonian. 



III. THE INTERACTION BETWEEN 
d-ELECTRONS 



This allows us to conclude that the nearest neighbor 
dxz-dxy and dyz-dxy mixing give rise to exz,xy{k) = 

-'2itxz,xysm{kx) and ^yz,xy(k) — i'ltyz^xySiniky). 

Now we can see how the 3-orbital model can reproduce 
the four-pocket Fermi surface. Let us begin by assum- 
ing that exz{k) lies just below the Fermi energy (-0.2eV) 
aX k = Y and disperses rapidly upward to wards M, 
reaching 2.5eV. From F, it descends towards X where its 
energy is — 1.4eV. The dispersion is relatively flat along 
X-V-Y and a shallow local maximum appear at F. The 
Fermi surface correspond to this band is shown in Fig. 
3a. The dyz band is similar except rotated by 90°. The 
reason for the choice of locating the dxz + dxy band at Y 
(as opposed to X) will be explained later. 

Next, we turn on the hybridization between dxz and 
dyz- The hybridization is maximal at X (midway be- 
tween X and Y) and creates two sheets which touches at 
F and two Fermi surfaces. We also find two concentric 
small hole pockets at F which become the well known 
hole pockets at F after folding. 

The dxy band ^xy{k) is assumed to be at — 0.5eV at 



Next, let us consider the single ion interaction between 
the electrons on the Fe d-orbitals which is given by 

Hi = \Y. I d'a^d^a;' ciix)c^{x)V{x - x')cl{x')c^,{x') 

a, a' 

Q-i ,0.2 ,0.3 ,0-4, a, a' 

d^xdi^x' (f)ai{x)(f)a^{x')V{x ~ x')(f>a3{x')(f)ai{x) 



where ai, • • • , 04 = xz, yz, xy label the orbitals, 0a (a;) is 
the wave function of the orbitals, and Cq is the electron 
operator for the a-orbital. Due to the symmetry of the 
orbitals, a,j must appear in pairs and the above can be 
rewritten as 



a a^b 
+ ^ JE ^aaclbCl3aCab + | E ^aacla'^UbCab (10) 



a^b 
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where Ui = Uaa = Ubb, U2 = Uab = Ui- 2J,^^ and 

j = J d'^xd'^x' Mx)Mx')v{x - x')Mx')Mx)- 

(11) 

Here the Ui term represents the Coulomb interaction 
between electrons on same orbital. The U2 term repre- 
sents the interaction between electrons on different or- 
bitals. The first J term is the exchange effect that favor 
parallel spins on the same Fe and is responsible for the 
Hund's rule. The second J term is pair hopping between 
different orbitals. 



IV. PAIRING INSTABILITY 

The above rf-electron interaction may cause a pairing 
instability. One possible pairing instability is the spin- 
triplet pairing between X and Y electron pockets.^'' In 
this paper, we will consider a different type of pairing 
- the spin-triplet p-wave pairing within the same Fermi 
pocket. 

First, let us consider the p-wave pairing on the Fermi 
pocket near the Y point (see Fig. 3). Such a Fermi 
pocket is a mixture of the dxz and dxy orbitals. Let ipxz 
and tpxy be the electron operators in the continuum limit 
in the dxz and dxy orbitals near the Y point. The electron 
operator ^ near the pocket at the Y point (see Fig. 2a) 
is a mixture of ipxz and Vxy 

The Hamiltonian has the Pr,Px, Pi-P,,, P^t,, PrT^., and 

-y 



PzTy symmetries, where Px x ^ —x and Py : y 
are reflections about a Fe atom. These symmetry dic- 
tates certain form of the continuum Hamiltonian. Alter- 
natively, we can expand the tight-binding picture near Y 
to obtain the following form which respect all the sym- 
metry requirements 

Ho = ihkl + hf4 + 4Mz^.^ (12) 
+ {iikl + hkfj + ely)^iy^Jxv + hkx{\^iyXl,xz + h.c.) 

Here k is the pseudo crustal momentum measured from 
the Y point. From Fig. 4(a), we see that the dxz band 
is flat along the y-axis, which implies S> ^2- Similarly 
^2 ^ ^1 • In the following, we will ignore <2 and ti . 
The electrons have two bands with energies 



(13) 



where 



£3 = ^{tikl - i2kl + - el 
€2 = takx- 



xyj 



xyj 



Ipxz and t/jxy are related to the electron operator ip in the 
upper band as 



ipxz{k) =Ukip{k), ipxy{k) =Vkip{k), (15) 



where 



Uk = 



Vk = 



162 



£3 - V4 + ^ 



^2el + 2el- 263^/4+^: 



(16) 



(14) 



To obtain the pairing interaction near the Y point in 
the spin-triplet channel we set a. = = a' = [3' =] in 
(10). We find that the first and the fourth terms in (10) 
vanish. The second and the third terms become 

^ Vy{k2Mm{k2)^l^^{-k2)]H^{kM^{-k^) (17) 



in the spin-triplet and (fe, — fe) pairing channel, where the 
effective pairing interaction of V' is 

V"y(fe2, fel) = -(J - U2)u*k^V*_k2UkiV-k, (18) 

From the effective pairing interaction, we can obtain a 
dimensionless coupling constant^* 

, Il^iIwR^l9^ik)Vik,k')gik') 

/(4?talP(^)P ^ ^ 

where J dak is the integration over the Y Fermi sur- 
face and Vk is the Fermi velocity. The function g{k) 
is a square harmonics which describes the shape of the 
superconducting gap, e.g. g{k) = 1 corresponds to an 
s-wave and g{k) = sin kx or sin ky a p-wave siipercon- 
ductor. The superconducting transition temperature Tg 
is given by Tc = Oe~^/^, where O is of order of the Fermi 
energy of the pocket: 17 ^ 0.2eV. 

From (18), we find that, when J > U2, the pairing 
interaction V^(fe2,fei) induces a pairing ^(fe) that have 
the same symmetry as MfcWfc. Because itfeffc is odd in kx 
(see Fig. 4b), the induced pairing is in p-wave channel 
g{k) = sinkx- Thus when J > U2, the rf-electron in- 
teraction will cause a spin-triplet p-wave pairing. Note 
that the location of the node depends on the choice of 
sinfca; versus sinfcj,, which in turn hinges on our assign- 
ment of the orbital to be X2;-like near Y . Triplet pairing 
has been proposed earlier, and fully gapped states such 
as Px + ipy were suggested. In contrast, after including 
the fe-dependent orbital-mixing, we find the on-site fer- 
romagnetic interaction to favor a particular nodal Px and 
Py states in the Y and X valleys respectively. 

As long as > e^z, we see from Fig. 4a that the Fermi 
surface changes its character from pure dxy to mostly dxz 
as a function of angles. |ufep and \vk\'^ must cross at 
some angles where |ufcWfc| takes its peak value 1/2 (see 



5 



Fig. 4b). Thus the effect of litfeiifel or A is relatively 
insensitive to doping provided that ep > e^^. Since the 
density of states is also independent of doping in 2D, this 
explains why is somewhat insensitive to doping and 
may extended to the hole doped side,^ except near zero 
doping. There the Ui term will drive an SDW instabil- 
ity due to the nesting between the electron and the hole 
pockets. ^^'-^^ 

We would like to mention that since E{k) = E{—k), 
our intra-pocket p-wave pairing appears even when the 
attraction J — U2 is weak. When the attraction J — 
U2 is strong enough to overcome the intcr-pocket energy 
splitting (~ 0.05eV), our model also has the instability 
in the spin-triplet inter-pocket pairing channel proposed 
in Ref. 17. Since the intra-pocket effective attraction 
is reduced by the matrix elements as shown in (18), the 
inter-pocket pairing may be stronger in large J—U2 limit 
while the intra-pocket pairing is stronger in small J — U2 
limit. 

Similarly, we can consider the spin-triplet pairing on 
the Fermi pocket near the T point (see Fig. 2b). Such 
a Fermi pocket is a mixture of the dxz and dyz orbitals. 
The dispersion and the d^z-dyz mixing near F can be 
determined from the symmetry consideration. The cru- 
cial difference is that the hybridization matrix element is 
now proportional to k^ky. After a similar calculation, we 
find that the spin-triplet pairing potential Vp(A;2, fci) sat- 
isfy Vf{k2,ki) = Vp(A;2,— fci) = Vf{—k2,ki). A p- wave 
pairing will result in a vanishing dimensionless coupling 
A = 0. Thus the Coulomb interaction in the rf-orbitals 
does not induce spin-triplet p-wave pairing on the two 
pockets near F even when J > 1/2- We see that the d^z- 
dxy (dyz-dxy) mixing at Y (X) in the 3-orbital model 
is crucial for the appearence of our p-wavc instability. 
A 2-orbital model has the wrong Fermi surface topol- 
ogy in that one hole pocket is located at F and M in Fig. 
2(b), instead of both being at F.^^ Although the 2-orbital 
model may allow certain superconducting states, it 
does not have the proper symmetry and the orbital mix- 
ing to generate the p-wave pairing proposed here. 



V. CONCLUSION 

In this paper, we examine the possibility that Coulomb 

interaction between the electrons in the Fc rf-orbitals may 
induce a superconducting phase. We find that when the 
Hund's rule ferromagnetic interaction on Fe is strong 
enough, ie when J > U2, the Coulomb interaction can in- 
duce pairing instability. In the weak coupling limit, the 
leading instability is found to be a spin-triplet p-wave 
pairing on the electron pockets at M (or at X and Y 
in the extended Brillouin zone). Although the Coulomb 
interaction can only directly induce pairing on the elec- 
tron pockets at M, the proximity effect will lead to a 
p-wave pairing on the hole pockets at F. So the resulting 
superconducting state has node lines on the 3D Fermi 
surfaces. 



The spin triplet pairing order paraneter is a compex 
vector d. The relative phases and relative orientations of 
the two order parameters and dy on the two pockets 
X and Y can have interesting relations which cannot be 
determined from the linear response calculation adopted 
here. One possible distribution of the phases is given 
in Fig. 2b, where p-wave pairing gap is positive in the 
shaded regions and negative in unshaded regions. 

At the atomic level, U « 3 to 4 eV and J « 0.7 eV, so 
U2 — J is positive. However, in a tight binding model in- 
volving only the Fe rf-orbitals, the As orbitals have been 
projected out and the appropriate U and J are those 
corresponding to the Wannier orbitals which are much 
more extended than the atomic d orbitals. A recent es- 
timate by Anisimov et al."^^ foimd the average U to be 
strongly renormalized down to 0.8 eV while J remains 
large at 0.5 eV. These are the more appropriate bare pa- 
rameters in a 3-band model. Furthermore, in a crystal, 
the strong hopping leads to extended quasiparticles near 
the Fermi surface. The on-site U2 and J will induce effec- 
tive interaction [/| and J* between those quasiparticles. 
The term induced by U2 will remain short ranged. How- 
ever, J may induce a long range couple because paral- 
lel spin configuration favors hopping, ie Hund's rule and 
hopping are compatible. Since the Fermi pockets are 
small, the effective interaction is given by ?7|(g) — J*{q) 
with q ^ kp ^/a- A long range J-coupling enhances 
J*{q = 0) and can potentially lead to a sign change and 
a net effective attraction. 

Since the original submission of this paper the exper- 
imental situation has evolved rapidly. Here we attempt 
a brief summary of the relevant experimental data. The 
issue of whether gap nodes exist remains open to debate. 
Photoemission data on Bai_,,;K2;Fe2As2 indicate an al- 
most isotropic gap on the hole pocket in this hole doped 
material. It was recently found that the magnetic field 
dependence of the specific heat in this material is linear, 
in contrast with the behavior taken as evidence for 
nodes in the electron doped material.^ At the same time, 
the nuclear spin relaxation rate ^ fits the law over 

three decades in La(Oi_a;Fex)FeAs.^^ Thus at the mo- 
ment, existing data seem to point to gap nodes in electron 
doped materials and their absence in hole doped materi- 
als. We note that, because the size of the Fermi pockets 
are small, it is usually advantageous to hide the nodes in 
the k space between Fermi pockets. In order to produce 
a gap node on the Fermi surface of the small pocket, one 
needs an effective pairing potential which varies rapidly 
on the scale of the small pocket. Our theory is one of the 
few that will do this, and we rely on the rapid k depen- 
dence of the hybridization matrix element. The message 
that wavefunction and matrix elements may play an im- 
portant role and must therefore be handled properly has 
validity beyond the special p-wave pairing scenario de- 
scribed here. 

What about singlet vs. triplet pairing? The Knight 

shift is probably the best way to answer this question. 
In a triplet superconductor, the spin contribution to the 
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Knight shift drops below Tc to zero if the magnetic field 
is parallel to the d vector, but remains unchanged if it 
is perpendicular. If the d vector is free to rotate, it will 
turn perpendicular to H and no change in Knight shift 
is predicted. This is apparently the case for Sr2Ru04. 
On the other hand, if d is locked to the latticx;, we ex- 
pect to see in a polycrystalline sample a drop of 2/3 of 
the value compared with singlet pairing. Without accu- 
rate knowledge of the orbital contribution to the Knight 
shift, this is hard to distinguish. Thus NMR on single 
crystals is needed to settle this question. As of this writ- 
ing, the only single crystal data available is from Ning 
et al?^ on BaFe1.8Coo.2As2. The data does not support 
triplet pairing in that a drop in the Knight shift is seen 
for field directions both parallel and perpendicular to the 



plane. On the other hand, a recent paper by Nakai et 
al?'^ on the FeP system Lao.srCao.isFePO shows that 
the magnetic behavior is quite different from the FeAs 
system. The Knight shift increases with decreasing tem- 
perature, indicative of ferromagnetic fluctuations above 
Tc and shows a very unusual increase below Tg. The 
authors speculate that magnetic fluctuations associated 
with triplet pairing may be responsible for the increase. 
The experimental situation remains in flux and it may 
be possible that different pairing scenarios may be com- 
peting and win out in different materials. 

This research is supported by DOE grant DE-FG02- 
03ER46076 (PAL) and by NSF Grant DMR-0706078 
(XGW). 
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